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High mobility two-dimensional hole system on hydrogen-terminated silicon 
(111) surfaces 
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We have realized a two-dimensional hole system (2DHS), in which the 2DHS is induced at an atomically flat 
hydrogen-terminated Si(lll) surface by a negative gate voltage applied across a vacuum cavity. Hole densities 
up to 7.5 x 10 11 cm -2 are obtained, and the peak hole mobility is about 10 4 cm 2 /Vs at 70 mK. The quantum 
Hall effect is observed. Shubnikov-de Haas oscillations show a beating pattern due to the spin-orbit effects, 
and the inferred zero-field spin splitting can be tuned by the gate voltage. 



Metal-oxide-semiconductor (MOS) field-effect tran- 
sistors (FETs) have long been used to study two- 
dimensional systems in silicon. While there has been 
extensive research on two-dimensional electron systems 
(2DESs) in silicon using MOSFETs over the last sev- 
eral decades, 1 ' 2 comparably little work has been done on 
two-dimensional hole systems (2DHSs), 2 5 primarily due 
to the fact that holes have lower mobility than electrons, 
and detailed study was not possible at such low mobility 
silicon hole devices. 

Recently, a new type of silicon FET device has been de- 
veloped, in which an atomically flat silicon (111) surface 
is terminated by a monolayer of hydrogen atoms, and a 
2DES is induced at the H-Si(lll) surface by a positive 
gate voltage through a vacuum barrier. 6 2DESs in these 
vacuum FET devices show extremely high mobility. 6 8 In 
this letter, we report that the vacuum FET technique can 
be extended to 2DHSs. A peak hole mobility of ~ 10 4 
cm 2 /Vs is obtained, exceeding by one order of magnitude 
the mobility in MOS structures and comparable with the 
mobility of pseudomorphic Si/SiGe heterostructures, 9,10 
where 2DHSs reside in the SiGe channels. The inte- 
ger quantum Hall effect (IQHE) is also observed. At 
lower magnetic fields, Shubnikov-de Haas (SdH) oscilla- 
tions show a beating pattern, which is due to the zero- 
field spin splitting. 11 In silicon, the spin-orbit effects lift 
the spin-degeneracy of the heavy hole subband even in 
a zero external magnetic field. The zero-field spin split- 
ting is from the structure inversion asymmetry (SI A), 12 
and can be tuned by the external gate. The zero-field 
spin splitting is of great interest because of its potential 
use in spintronics devices 13 and in studying fundamental 
physics. 12 

The vacuum FET device includes two pieces, as shown 
in Fig. 1. One is a p~ Si(lll) piece (float zone, resistance 
> 10000 (]-cm), and the other is a silicon-on-insulator 
(SOI) piece. First, four p + contact regions are formed 
in the Si(lll) piece by ion implantation with 2.4 x 10 15 



2DHS onH-Si(111) 



p- Si(111) 



(drain) 



SOI- 



(shie 
■ I P-Si(K 



(shield) 

p+ Sf (gate)' 
Si(100) 




a ) To whom correspondence should be addressed. Electronic mail: 
hubh@lps.umd.edu 

b ) Present address: LSW, Arlington, VA 22203, USA. 



FIG. 1. (Color online) Schematic cross section of a vacuum 
FET (not to scale). A H-Si(lll) piece is contact bonded to 
a SOI piece in a vacuum. The H-Si(lll) piece has heavily 
boron-doped (p + ) regions to make contact to the 2DHS. In 
the SOI piece, the shield and the gate are p + layers formed by 
double boron ion implantation. Dry etching is used to create 
the cavity. The arrows depict the electric field. The 2DHS is 
induced at the H-Si(lll) surface in the encapsulated vacuum 
cavity. 



cm -2 , 15 keV boron ions through a 30 nm thick thermal 
oxide, as shown in the insert of Fig. 2. It is annealed in 
N2 gas at 1000° C for 10 minutes to activate dopants and 
reduce defects. Second, in the SOI piece, the shield and 
gate conducting layers are formed by double boron ion 
implantation, and a cloverleaf-shaped cavity is created 
by dry etching. Finally, the Si(lll) piece is H-terminated 
by immersing it in an ammonium fluoride solution, and 
then these two pieces are contact bonded in a vacuum 
chamber. The detailed fabrication processes and the ba- 
sic operating principle have been discussed elsewhere. 6,14 
Here the SOI piece acts as a remote gate to induce a 
2DHS at the Si(lll) surface through the encapsulated 
vacuum cavity. This cavity also protects the air sensitive 
H-Si(lll) surface. Compared to the Si/Si02 interface in 
MOSFETs, the vacuum/H-Si(lll) interface is atomically 
flat and has much fewer defects, leading to a much higher 
carrier mobility. 

The devices are characterized in a dilution refrigera- 
tor with a base temperature of ~ 70 mK. Longitudinal 
(R xx ) and Hall (R X y) resistances are determined by Van 
der Pauw measurements with the standard low-frequency 
AC lock-in technique using a 7 Hz, 100 nA current source. 
As shown in the insert of Fig. 2, the Van der Pauw geom- 
etry consists of a center cloverleaf-shaped 2DHS, which is 
electric field induced, and four heavily boron-doped con- 
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FIG. 2. (Color online) Left axis: mobility as a function of gate 
voltage at 70 mK and 4.2 K. At 70 mK, the wild fluctuations 
of the mobility are due to the high contact resistances to the 
2DHS. Right axis: Hole density vs gate voltage. The hole 
densities (triangles) are determined by the magnetic field B 
at the filling factor v = 5 minima: p2d = veB/h. The dashed 
line is a linear fit. Insert: on the Si(lll) piece, Van der Pauw 
geometry consists of a cloverleaf-shaped 2DHS with a 0.5 x 0.5 
mm 2 center square, and four heavily boron-doped contacts. 



tacts (labeled as 1, 2, 3, 4). Four-terminal resistance is 
denned as it^,j m = Vimlhj (ij,l,m = 1,2,3,4), where 
current Iij is fed from contact i to contact j, and volt- 
age Vim is measured between contacts 1 and m. The 
sheet resistance is determined by the standard Van der 
Pauw technique. 15 Although the Van der Pauw geometry 
is symmetric, R xx and R yy are found to differ by ~ 35%, 
so the sheet resistance, as well as the mobility, is the aver- 
age of x and y directions. 16 The anisotropy can be caused 
by some mechanisms such as miscut, disorder, strain and 
the non-parabolicity of the valence band, 7,8,17 but the 
exact sources are yet to be ascertained. Hole density p2d 
is determined by the magnetic field B at the filling factor 
v = 5 minimum of the longitudinal resistance curve at T 
~ 1 K: p2d = veBjh^ where e is the electron charge and 
h is Planck's constant. The data (triangles) are shown in 
Fig. 2, and the dashed line is a linear fit. Hole densities 
up to 7.5 x 10 11 cm -2 (V g = -110 V) are obtained. If 
we use a parallel plate capacitor model, the equivalent 
depth of the cavity is 759 nm, which is consistent with 
the measurement result from a profilometer. The extrap- 
olated threshold voltage is —3.9 V, which indicates the 
density of the trapped charge at the surface ~ 2.8 x 10 10 
cm -2 . The mobility of the 2DHS is also shown in Fig. 
2, and the peak mobility is about 9800 cm 2 /Vs at gate 
voltage V g = —100 V, T = 70 mK on this device. This is 
one order of magnitude higher than the hole mobility of a 
Si(lll) MOSFET. 4 The contact resistances between the 
contacts and the 2DHS increase dramatically at \V g \ < 60 
V, when the temperature is lowered from 4.2 K to 70 mK. 
For example, they increase from ~ 20 KQ to ~ 400 KQ 
at V g = —50 V. It is difficult to get reliable measure- 
ments at this region, and the calculated mobility shows 
spurious wild fluctuations at T = 70 mK. Our discussion 



FIG. 3. (Color online) Longitudinal resistances i?i3,24, i?i2,34, 
and Hall resistance R xy vs perpendicular magnetic field B at 
V g = -90 V (p 2d = 6.27 x 10 11 cm" 2 ) and T = 70 mK. 
Quantum Hall effect is observed near B = 5.2 T. The insert 
shows an enlarged section of i?i3,24, which exhibits a beating 
pattern, and the arrow marks the node location (B = 1.14 
T). 



below will focus on the data at \V g \ > 60 V. 

As evidence of the high mobility 2DHS, Fig. 3 shows 
the longitudinal (R xx ) and Hall (R X y) resistances as a 
function of perpendicular magnetic field B at V g = — 90 
V (p 2 d = 6.27 x 10 11 cm -2 ) and T = 70 mK. The integer 
quantum Hall effect (IQHE) is observed. At the vicinity 
of B = 5.2 T, R xx approaches zero, and R xy develops 
a plateau at v = 5. Clear SdH oscillations appear at 
B field above 0.4 T. A close look at the trace of R xx 
reveals a beating pattern of the SdH oscillations with a 
node near B = 1.2 T, as shown in the insert of Fig. 3. 
Although the beating pattern has been observed in p- 
channel Si(110) devices, 3 ' 5 to our best knowledge, this is 
the first time it has been observed on a Si(lll) surface. 
The beating pattern can be explained by the spin-orbit 
effects, although different explanations exist. It could be 
from gross inhomogeneities of the carrier density in the 
sample. 18 Two samples from different wafers have been 
measured, and both show a similar beating pattern. At 
V g = -90 V, the beating node is located at B = 1.14 T 
for the device reported here, and it is at B = 1.17 T for 
the other device. It is therefore highly unlikely that the 
gross inhomogeneities are the cause. 

The beating pattern indicates that there are two 
closely spaced frequency components with similar am- 
plitudes, arising from two populations of holes. The 
hole subband densities p\ d and p^ d can be determined 
by the Fourier analysis of the SdH oscillations with 
p 2d = g s ef ± /h, where g s is the degeneracy factor, /+ 
and /~ are the frequencies of the SdH oscillations (Fig. 
4(b)). 19 ' 20 Here R xx vs B~ x data in the range 0.2 < B < 
4 T are Fourier transformed after subtraction of a slowly 
varying background, and peak locations are used as the 
frequencies of the SdH oscillations. For \V g \ > 85 V 
(with two well-resolved peaks), g s must be 1 in order 
to have the sum of the subband densities equal to the 
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FIG. 4. (Color online) (a) Ai?i 3 ,2 4 (= #13,24 (5) - #13,24 (0)), 
A#i2,34(= #12,34(5) — #12,34(0)) vs filling factor v at Vg = 
—90 V and T = 70 mK. The corresponding magnetic field B 
is shown on the top axis. The node location (B = 1.14 T) 
is marked by an arrow. Across the node, the minima change 
from the odd numbers of v to the even numbers of v. (b) 
Fourier power spectra of the SdH oscillations (#13,24) at dif- 
ferent gate voltages. Maximum peak amplitudes are normal- 
ized to 1. At V g — —90 V, the solid line is the power spectrum 
of #13,24, and the dotted line is the power spectrum of #12,34- 
(c) Zero-field spin splitting and total spin splitting (both di- 
vided by eh/m*) vs gate voltage. Zero- field spin splitting is 
determined by Fourier power spectrum (b). The total spin 
splitting is determined by the node location (a), and the er- 
ror bars are estimated by the error bound [is no de — 1, Vnode + 1] , 
where v no de is the node location in terms of filling factor. 



total density, so these two components can only be from 
one subband with two spin-split levels, i.e. two spin-split 
subbands. 

Landau levels of each spin population give rise to 
magnetoresistance oscillations with similar amplitude, 
described by cos[27r(£^/ ± S/2)/huo c ], where Ef is the 
Fermi energy, S is the energy separation between the 
two spin-split subbands, H is ft/27r, and uo c is the cy- 



clotron frequency. 11 The sum of these two components 
leads to the modulation of the SdH oscillations given 
by A ~ cos(7t5/Huj c ). Beating nodes are located at 
5/huo c = n = ±0.5, ±1.5, etc. Only one node is clearly 
identified near B ~ 1.2 T (Fig. 3). Fig. 4(a) shows the 
SdH oscillations in terms of filling factor v at V g = — 90 
V. At v < 22.7, the minima of the SdH oscillations are 
at the odd numbers of v. At v > 22.7, the minima are 
at the even numbers of v. This manifests the transition 
across the node. We assume that the node corresponds 
to n = 0.5 (highest-field), because the n = 0.5 node 
is most easily observed, and there is no other node up 
to 6 T. 13 The n = 1.5 node is estimated to be around 
B ~ 0.4 T, one third of the n = 0.5 node, so it is dif- 
ficult to observe for this device. From the node loca- 
tion B no de, the total spin splitting can be determined by 
S = 0.5huj c = (eh/m*)(0.5B no( ie), where m* is the effec- 
tive hole mass. In Fig. 4(c), the total spin splitting S is 
described by m*5/eh = 0.5B no d e - 

In theories of the spin-orbit interaction, spin split- 
ting can have complicated B dependence. 12 The total 
spin splitting S can be expanded as 8 = So + Si(Huj c ) + 
52(huj c ) 2 + where So is the zero-field spin splitting, 
S\{hjj c ) is the linear splitting, etc. 11 For low magnetic 
fields, if only the first two terms are considered, then 
A ~ cos(7r£o/ftco> c + ttSi). The frequencies of the SdH 
oscillations can be used to determined the zero-field spin 
splitting by So = (efi/m*)(/ + — /~). So the zero-field 
spin splitting So can be described by m*So/eh = / + — / _ , 
also shown in Fig. 4(c). There are two important features 
in the data. First, the zero-field spin splitting increases 
with increasing perpendicular electric field E z . Second, 
at | V g \ > 85 V, So > S, which means that a perpendic- 
ular magnetic field can reduce the spin splitting. More 
detailed analysis is needed to better appreciate the re- 
sults. These topics will be explored in future studies. 

In conclusion, the vacuum FET technique has been 
demonstrated to create high mobility 2DHSs in Si(lll). 
In these high mobility 2DHSs, we are able to observe the 
IQHE, and the beating pattern of the SdH oscillations on 
Si(lll) surfaces, which are difficult to reach previously. 
There are many interesting questions yet to be resolved, 
such as the non-parabolic valence band structure and the 
spin-orbit effects. Our high mobility 2DHSs open up new 
opportunities to explore fundamental physics in the two- 
dimensional hole systems. In addition, with the ability 
to create both a 2DES and a 2DHS on a Si(lll) surface, 
it is straightforward to develop a bipolar surface device. 
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